Example 7.9

In [2]:

import numpy
import matplotlib.pyplot as plt
%matplotlib inline

In [16]:

omega = numpy.logspace(-3, 1, 1000)
s = 1j*omega

In [17]:

Kc = 1.13
G = 3*(-2*s + 1)/((5*s + 1)*(10*s + 1))
def K(Kc):
    return Kc*(12.7*s + 1)/(12.7*s)
Gprime = 4*(-3*s + 1)/(4*s + 1)**2
w_I = (10*s + 0.33)/((10/5.25)*s + 1)

In [18]:

K1 = K(Kc=1.13)
T1 = G*K1/(1 + G*K1)

In [19]:

K2 = K(Kc=0.3)
T2 = G*K2/(1 + G*K2)

In [20]:

plt.loglog(omega, numpy.abs(T1), '-', label='$T_1$ (not RS)')
plt.loglog(omega, numpy.abs(T2), '-', label='$T_2$ (RS)')
plt.loglog(omega, numpy.abs(1/w_I), '--', label='$1/W_I$')
plt.legend()

Out[20]:

<matplotlib.legend.Legend at 0x10d7bcdd0>
../_images/notebooks_Example_7.9_6_1.png 
In [ ]:

e